1/ / noise in experimental Sinai quantum billiards 
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It was recently conjectured that 1 // noise is a fundamental characteristic of spectral fluctuations 
in chaotic quantum systems. In this Letter we show that the level fluctuations of experimental 
realizations of the Sinai billiard exhibit 1/f noise, corroborating the conjecture. Assuming that 
the statistical properties of these systems are those of the Gaussian orthogonal ensemble (GOE), 
we compare the experimental results with the universal behavior predicted in the random matrix 
framework, and find an excellent agreement. The deviations from this behavior observed at low 
frequencies can be easily explained with the semiclassical periodic orbit theory. In conclusion, it is 
shown that the main features of the conjecture are affected neither by non-universal properties due 
to the underlying classical dynamics nor by the uncertainties of the experimental setup. 
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Quantum chaos is a well established subject since the pioneering work of Bohigas and collaborators|l| (for further 
review see, for example, 0,0). Recently, a different approach to quantum chaos has been proposed jjT Considering 
the sequence of energy levels as a discrete time series, and the appropriate statistic related to the fluctuations of 
the excitation energy, it was conjectured that chaotic quantum systems are characterized by 1/f noise, whereas 
integrable ones exhibit 1/f 2 noise. Since chaotic spectra behave like antipersistent time series, and regular ones as 
neither persistent nor antipersistent time series, it was possible to establish a relationship between spectral rigidity 
and antipersistence. This conjecture is supported by numerical calculations which involve atomic nuclei and the 
classical Random Matrix Ensembles (RME) Q. Later, theoretical expressions were derived in Ref. [f| using random 
matrix theory. They explain, without free parameters, the universal behavior of the excitation energy fluctuations 
power spectrum. The theory gives excellent agreement with numerical calculations. It also reproduces the 1/f power 
law characteristic of chaotic systems and the 1/f 2 law characteristic of integrable system through several orders of 
magnitude of the frequency interval. 

All these results were obtained in the Random Matrix Theory (RMT) and Nuclear Shell Model frameworks. These 
models allow an accurate statistical analysis since very long sequences can be obtained with high precision; moreover, 
RMT is the usual theoretical model in studying quantum chaos, and the high excited states of the atomic nuclei follow 
RMT predictions almost perfectly^ Q. However, other quantum chaotic systems cannot provide such quality of data. 
Experiments usually provide short level sequences with some precision problems: if two resonances are separated by 
a distance smaller than the experimental line width, they are registered as one; in microwave billiard experiments, 
certain resonances are missing if the emitting antenna is positioned close to a nodal line. Quantum systems with 
a classical analogue can provide very long sequences, but they present some non-universal features that spoil the 
universal behavior of long-range correlations ja, 0J , which reveal some dynamic properties of the classical analogue 
-the behavior of the shorter periodic orbits of the system. Therefore, analyzing whether experimental limitations and 
the existence of short periodic orbits modify our previous result becomes an important issue. 

In this Letter we study the spectral fluctuations of experimentally modeled quantum Sinai billiards. We find that 
the excitation energy level fluctuations of these systems are in clear agreement with the theoretical expressions found 
in Ref. Therefore, they exhibit 1/f noise, in spite of small deviations due to experimental limitations and non- 
universal properties due to short periodic orbits. Moreover, we present a simple analysis of these anomalies using 
periodic orbit theory. 

In the present work, quantum billiards have been substituted by suitably shaped microwave resonators, whose 
eigenfrequencies are measured. This kind of experiment uses the fact that the time-independent Schrodinger equation 
and the wave equation are mathematically equivalent. The data we analyze in this work has already been used in 
Refs. [ToL ITT| . although in a different context. 

To be specific, we have used a quarter Sinai billiard varying the values of the length I as depicted in Fig. ^ The 
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FIG. 1: Sketch of the quarter Sinai billiard with varying length printed in scale. The length of the billiard was varied from I 
= 480 mm to 460 mm in steps of 0.2 mm. The width b — 200 mm and the radius of the quarter circle r = 70 mm. 



experimental setup makes it possible to measure resonances in the frequency range 14.5 to 15.5 Ghz, where there 
are about 120 eigenvalues for each value of the parameter I. We have considered about 100 of such sequences as 
realizations of a stochastic process representative of the billiard, as we have previously done in numerical experiments 

B- . 

Prior to the statistical analysis, we have to remove the modulation induced by the mean level density in the spectra. 
This procedure is called unfolding; it consists in locally mapping the actual energy levels Ei into new dimensionless 
levels £j 

E % — >e, =N(Ek), iE [l,--- ,N], (1) 

where N is the dimension of the spectrum and N(E) is a smooth approximation to the actual accumulated density. 
Quite generally, N(E) can be calculated with the aid of the following rule: in the limit of larg e E, the average density 
of eigenstates is proportional to the volume of the phase space with energy less than E [13 . For quantum billiards 
we can use an improved version of this result, which is also valid near of the ground-state [9j: 

N(E) = N Weyl {E) + N bb (E) . (2) 

Weyl 

N is the so called Weyl law; it includes a main term related to the volume of the phase space, and two correction 

terms related to the surface and the curvature of the boundary . N (E) is a further correction term which takes 
into account the smooth trend produced by a set of non-isolated marginally stable periodic orbits bouncing between 
the two straight sides of the billiard. 

After performing the unfolding procedure, we can calculate the nearest neighbor level spacings defined by Si = 
ti+i — ei, i G [1, • • • , N — 1]; since the unfolded spectrum has a mean level density equal to 1, their average value is 
(s) = 1. To study the long-range correlation structure, we use the 5 n statistic 0,03 defined by 

n 

5 n = y~\(si ~ (»» = £n+i - £i - n, ne [0, ••• ,N- 1]. (3) 



This statistic represents the deviation of the excitation energy of the (n+l)ih unfolded level from its mean value n, 
and it is closely related to the spectrum fluctuations. Indeed, we can write 5 n = —N(E n +i) if we appropriately shift 
the ground state energy; thus, 6 n represents the cumulated level density fluctuations. If we consider the index n as the 
analogue of a discrete time |l6j , we can analyze S n with numerical techniques normally used in the study of complex 
systems. The most simple is the analysis of the power spectrum obtained from the discrete Fourier transform (DFT) 
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S (iv k ) = 5 k , (5) 
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where 

27rfc , . 

Uk = NA (6) 

and A is the sampling interval of the signal. 

Here we face the following problem. The experiment cannot provide all energy levels in the interval 14.5 and 15.5 
GHz: there are a few missing eigenvalues due to nodal lines of the wavefunction at the position of the antenna. 
Nevertheless, due to the level dynamics measurement, one can detect where a certain level is missing. A compar- 
ison with eq. |J2J shows that all levels are accounted for. Thus, we are able to arrange the known energy levels 
as {e Pl ,e P2 ,- • • ,e Pi ,- ■ ■ ,£ PM }, where M{< N) represents the total number of levels whose energy is known, and 
{pi,P2, • • • ,Pm} are their positions in the spectrum. Therefore, the 8 n function can only be calculated for certain 
whole numbers. Taking e Pl as a fictitious ground state, we have 

Sn 3 = £ Pj+1 - £ P1 - rij , rij = pj+i - pi, j G [0, • • • , M - 1], (7) 

Consequently, S nj is not a standard time series, but an unevenly spaced time series. 

There are some obvious ways to get an evenly sampled signal from the unevenly spaced pj points. Interpolation 
is one possible way, but generally speaking, interpolation techniques perform poorly. Therefore, we have decided 
to proceeded in a different way to analyze the experimental data. For an arbitrarily sampled time series the power 

spectrum can be estimated by the periodogram, S(uj) = S(lo) , where S(u>) is given by 



S(u)) = S nj exp (—iuirij) . (8) 
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Eq. JSJ) is defined for any frequency ui. However, since the original signal is discrete, one should be able to chose a 
finite set of frequencies to evaluate the Fourier transform, and then using the values of eq. I© at these frequencies, 
recover the original signal. For evenly spaced time series of length N eq. © provides a natural set of frequency values. 
This set has been built from a fundamental frequency loi — 2n/(NA), which defines a sine wave of period equal to 
the whole interval length; the remaining frequencies are multiples of it. As commented above, the set {5(LOk)}, where 
fee [0, 1, • ■ • , N — 1], is known as the Discrete Fourier transform of the signal. 

In our case we have generated an appropriate set of frequencies as follows: 

(i) Since we know the smooth part of the level density (even if there are some missing levels) , we can perform the 
unfolding of the spectrum, leading to a constant average level density. Thus, we can set the fundamental sampling 
interval to A = 1. 

(ii) Given a window j\f < N, we select all the experimental billiards for which there is at least a pair of measured 
levels at positions p a and pb, satisfying that j\f — pt — p a - 

(iii) We consider the following frequency set uok = 2nk/Af with k = 1,2, This is probably the best choice 
because the fundamental frequency, io\ = 2ir /j\f, corresponds to a sine wave of period equal to the whole energy 
interval e Pb — e Pa ~ j\f ~ pb — p a - (Remember that after the unfolding procedure < s >= 1.) 

Due to the definition of our periodogram and to the selected frequency set, we have that S(ui + 2ir) = S(ui — 2ir) = 
S(lo). For this reason the power spectrum is defined in the basic interval [0, n). For evenly sampled series the frequency 
lu = 7r is known as the Nyquist frequency and a change of trend in the power spectrum can be observed near this 
upper limit. 

(iv) Finally, the average power spectrum of the S n signal is calculated using the selected billiards, and the experi- 
mental results are compared with the theoretical expression found in Ref . Q . Assuming that j\f ^ 1 , and that the 
statistical properties of the spectral fluctuations are those of the Gaussian orthogonal ensemble (GOE), the standard 
model for time-reversal invariant systems, the universal behavior of the S n power spectrum is given by 



(SH)) - — ^ — + ; 2 y + — -j^ + e, c fe e (o, tt), (9) 



Here. 



K M I 2r-rlog(2r + l),r<l 

K GO e(t) = \ 2 _ t log( |i±i )) T > ! (10) 
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is the spectral form factor of GOE, whose analytical expression is well known |14| . The constant shift O reflects that 
6 n is not a continuous two-point measure (defined for every energy interval) , but a discrete function (defined only for 
each energy level). For GOE O = —1/12, provided that the whole sequence of S n values is included in the calculation 
of the power spectrum. In this case, we cannot evaluate S n for all n, but only for some energy levels. This feature 
entails that the value of should be affected by the number M of levels whose energy is known. When uj k -C n, eq. 
© reduces to 

(S(u k )) = — , LU k « 7T, (11) 
TTWfc 

showing that for small frequencies the excitation energy fluctuations exhibit 1// noise in GOE-like systems. 
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FIG. 2: Averaged power spectra of the 5 n function for the Sinai microwave billiards (filled circles) compared with the 
theoretical prediction (solid line) as given by eq. The dashed straight line is the 1// noise predicted by the approximated 
expression (I I U . Vertical arrows indicate the frequencies of the two shortest periodic orbits of the billiard. 



In the forthcoming analysis we have chosen J\f = 85 since it is large enough to obtain significant results, and it 
is possible to find measured levels at positions p a and pb fulfilling pi, — p a = 85 in almost every sequence of levels; 
in fact there are 97 sequences satisfying this condition. Using a log-log scale, Fig. [2] compares the average power 
spectrum of 5 n and the theoretical predictions of eq. JHJ with w —1/38. This value was obtained by means of a 
least-squares fit to the numerical data. In spite of the experimental limitations, the agreement between the numerical 
and experimental results is very good. Moreover, it is clearly seen that the spectral fluctuations of these systems show 
a 1/ f power law. The significant power damping observed at lower frequencies can be understood by a simple analysis 
of the contribution of the bouncing ball orbits to the power spectrum. If we consider < I > = 470 mm and < b > = 
200 mm as the average length of the straight sides of the billiards, the frequencies corresponding to these orbits are 
u> = 0.09 and to = 0.21. These values define approximately the interval where the power is significantly diminished; 
therefore, we can conclude that this phenomenon is due to the removal of the contribution of bouncing ball orbits to 
the level density. All this is consistent with the idea that long period classical orbits give rise to an universal behavior 
for integrable and ergodic systems regardless of the system specificities |18| . Only when the main contribution comes 
from the scarce short periodic orbits, deviations from the universal behavior might be observed. 

In summary, we have studied the spectral fluctuations of quantum Sinai billiards modeled experimentally . We 
have considered the S n statistic, that characterizes the fluctuations of the excitation energy, as a discrete time series. 
Using about one hundred different lengths of the straight sides of the billiard, the average power spectrum of S n was 
calculated. Recall that performing an average is very important to reduce the fluctuations and clarify the secular 
trend of the 6 n function. The main result of this analysis is that the spectral fluctuations of these systems exhibit a 
clear 1/ / noise through almost the whole frequency domain. We have also compared the theoretical predictions of 
eq. Q for GOE-like systems with these results, finding an excellent agreement except in the lower frequency region. 
Nevertheless, a simple study reveals that these deviations from the universal behavior predicted by RMT are due 
to the removal of the bouncing ball orbits during the unfolding procedure; therefore there is no contribution to the 
power spectrum in the low frequency domain. Last but not least, we can say that S n statistic is not affected by the 
experimental limitations, therefore we can conclude that it is a quite strong statistic. 
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Prior to this work, we have found 1// noise in energy spectrum fluctuations for random matrix ensembles and 
atomic nuclei; these results allowed us to conjecture that 1// characterize spectral fluctuations of chaotic quantum 
systems. The new results obtained in this Letter corroborate this conjecture and give us a hint on the non-universal 
features that can appear in the power spectrum at lower frequencies due to the short periodic orbits. 
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